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A magnetic nanoparticle in a vortex state is a promising candidate for the information storage. 
One bit of information corresponds to the upward or downward magnetization of the vortex core 
(vortex polarity) . Generic properties of the vortex polarity switching are insensitive of the way how 
the vortex dynamics was excited: by an AC magnetic field, or by an electrical current. We study 
theoretically the switching process and describe in detail its mechanism, which involves the creation 
and annihilation of an intermediate vortex-antivortex pair. 
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Introduction 

Artificial mcsoscopic magnetic structures provide now 
a wide testing area for concepts of nanomagnetism and 
numerous prospective applications 0, 0. Investigations 
of magnetic nanostructures include studies of magnetic 
nanodots, which are submicron disk-shaped particles, 
which have a single vortex in the ground state due to 
the competition between exchange and magnetic dipolc- 
dipole interaction pj]. A vortex state is obtained in nan- 
odots that are larger than a single domain whose size is a 
few nanometers: e.g. for the Permalloy (NigoFe2o) nan- 
odot the exchange length I f« 5 nm. Having nontrivial 
topological structure on the scale of a nanomagnet, the 
magnetic vortex is a promising candidate for the high 
density magnetic storage and high speed magnetic ran- 
dom access memory For this one needs to control 
magnetization reversal, a process in which vortices play 
a big role Q . Great progress has been made recently with 
the possibility to observe high frequency dynamical prop- 
erties of the vortex state in magnetic dots by Brillouin 
light scattering of spin waves [|, time-resolved Kerr 
microscopy 0, phase sensitive Fourier transformation 
technique ff X-ray imaging technique 0], and micro- 
SQUID technique JhJ. 

The control of magnetic nonlinear structures using an 
electrical current is of special interest for applications 
in spintronics jllL The spin torque effect, which is 
the change of magnetization due to the interaction with 
an electrical current, was predicted by Slonczewski ]T^ ] 
and Berger [Q in 1996. During the last decade this ef- 
fect was tested in different magnetic systems (u| |l^] 
and nowadays it plays an important role in spintronics 
0, |||. Recently the spin torque effect was observed in 
vortex state nanoparticlcs. In particular, circular vortex 
motion can be excited by an AC or a DC |^0| spin- 
polarized current. Very recently it was predicted theoret- 
ically |2l], ^2) and observed experimentally |2j| that the 
vortex polarity can be controlled using a spin-polarized 
current . This opens up the possibility of realizing electri- 



cally controlled magnetic devices, changing the direction 
of modern spintronics [ p4[ . 

The basics of the magnetic vortex statics and dy- 
namics in Hciscnberg magnets was studied in 1980s, in 
particular by the group of Kosevich 
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• Typi- 
cally, the vortex is considered as a rigid particle without 
internal degrees of freedom; such a vortex undergoes a so- 
called gyroscopic dynamics in the framework of Thiele- 
like equations 
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|32j ] . Using a rigid vortex dynamics 
a number of dynamical effects were studied, for a review 
see Ref. [[33) . Cycloidal oscillations around the mean vor- 
tex trajectory can be explained by taking into account 
changes of the vortex shape due to its velocity [||[ Q . 

The rigid approach also fails when considering the vor- 
tex dynamics under the influence of a strong or fast exter- 
nal force. In particular, it is known that external pump- 
ing excites internal modes in vortex dynamics in Hciscn- 
berg magnets, whose role is important for understanding 
the switching phenomena [3(| |3^. [38| . |3fjfl , and the 
limit cycles in vortex dynamics JlGj, |4l[ . 

Very recently we have found that under the action of 
a magnetic field, or, alternatively, under the action of a 
spin current on a vortex state nanodot, the vortex core 
magnetization (polarity) can be switched. This effect 
was studied in Ref. [pi) analytically and confirmed nu- 
merically by direct spin-lattice simulations for Heisen- 
berg magnets. However, the presence of the long-range 
dipolar interaction crucially changes the vortex dynamics 
p2[ , because this interaction creates an effective nonho- 
mogeneous anisotropy ^l). The goal of this paper is to 
consider the generic properties of polarity switching in 
magnetic nanodots. 



MODEL AND EQUATIONS OF MOTION 

The magnetic energy of nanodots consists of two parts: 
Heisenberg exchange and dipolar interactions 
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Here S n is a unit vector which determines the spin 
direction at the lattice point n, £ = y/ A/ (fioMg) is 
the exchange length (A is the exchange constant, /xq is 
the vacuum permeability, Ms is the saturation magne- 
tization), the vector S connects nearest neighbors, and 
e nn > = (n — n')/\n — n'\ is a unit vector. The lattice 
constant is chosen unity length. 

It is known that as a result of the competition between 
the exchange and the dipolar interactions the ground 
state of a thin magnetically soft nanodisk is a vortex 
state: the magnetization lies in the XY disk plane in the 
main part of the disk and is parallel to the disk edge, 
forming the magnetic flux-closure pattern characterized 
by the vorticity q = +1. At the disk center the magneti- 
zation distribution forms a vortex core, which is oriented 
either parallel or antiparallel to the z-axis. The former 
is characterized by a polarity p = +1 and the latter by 
p = -l. 

The spin dynamics of the system is described by the 
modified Landau-Lifshitz-Gilbert equation 

— aS n x S n + S n x b n . (2) 
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Here the overdot indicates the derivative with respect to 
the dimcnsionlcss time r = u^t with ujq = ^wyMs, a<l 
is a damping coefficient. The last term is a spin torque 
due to external forces, which results in a magnetic field 
b n . In the presence of the homogeneous rotating field 
B{t) = (B cos lot, B sinwr, 0), the dimensionless field 
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Here and below all frequencies are measured in units of 
wo, and all distances in units of t. 

When an electrical current is injected in the pillar 
structure, perpendicular to the nanodisk plane, it influ- 
ences locally the spin S n of the lattice through the spin 
torque T n = S n x b n |S| Q, where 
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Here j = J e /J p is a normalized spin current, J e 
is the electrical current density, J p 
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is the disk thickness, 
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electron charge, 
3 = (1 + 

<E (0; 1) denotes 
the degree of the spin polarization; a = ±1 gives two 
directions of the spin-current polarization. 
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VORTEX STRUCTURE AND RIGID VORTEX 
DYNAMICS 

In the case of weak dipolar interactions, the char- 
acteristic exchange length £ is larger than the lattice 
constant a, so that in the lowest approximation in the 
small parameter ajl and weak gradients of magnetiza- 
tion m = —(Sn) = (sin cos (j>, sin 9 smcj), cos 9) we can 
use the continuum approximation for the Hamiltonian 
(0). 
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where H' 



is the normalized magnetostatic field, which 
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comes from the dipolar interaction 42 . In (|4|) we suppose 
that the magnetization distribution does not depend on 
the z-component (which is valid for thin disks) and we 
have normalized the energy by the disk thickness h. 

Let us consider a cylindric nanoparticle of top surface 
radius L and thickness h. For small size nanoparticles the 
ground state is uniform; it depends on the particle aspect 
ratio e = h/(2L): thin nanodisks are magnetized in the 
plane (when e < e c « 0.906 (4^]) and thick ones along 
the axis (when e > e c ). When the particle size exceeds 
some critical value, typically (3 — A)i [Q, the magne- 
tization curling becomes energetically preferable due to 
the competition between the exchange and dipolar in- 
teractions. For a disk-shaped particle there appears the 
vortex state. For thin enough disks the vortex structure 
does not depend on the z-coordinate of the disk and the 
magnetization distribution for the vortex, which is situ- 
ated in the disk center, is determined by the expressions 
11 

cos0=p/(|C|), = 9 argC + CV2. (5a) 

Here C, = x+iy is the coordinate in the disk plane, q = +1 
is the vortex -n\ topological charge (vorticity), C = ±1 
characterizes the vortex chirality, and p = ±1 determines 
the direction of the vortex core magnetization (polarity) . 
The bell-shaped function /(|CI) describes the vortex core 
magnetization. The vortex polarity is connected to the 
7T2 topological properties of the system, the Pontryagin 
index 

Q = — J c^xEijm ■ [djin x djin] . 

For the vortex configuration the Pontryagin index takes 
the half-integer values Q = —pq/2. 

For the vortex solution localized at Z = X + iY = 
i?exp(«4>), the (/>-field has the form </> = qa,rg(( — Z) + 
CV/2. Such a form of the </>-field satisfies the Laplace 
equation, so it describes the planar vortex dynamics (ne- 
glecting the z-component of magnetization) for the in- 
finite Hciscnberg magnets. The finite size effects for 
the circular Heisenberg magnet can be described by the 
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image- vortex ansatz jB3|]. Moreover, one has to take into 
account the long-range dipolar interactions, which lead 
in general to integro-differential equations fig} ]. For a 
thin ferromagnet the dipolar interactions produce an ef- 
fective uniaxial anisotropy of the easy-plane type caused 
by the faces surface magnetostatic charges |47| and 
an effective in-plane anisotropy caused by the edge sur- 
face charges (surface anisotropy) |l7], 48 1. Due to the 
surface anisotropy the magnetization near the disk edge 
is constrained to be tangential to the boundary, which 
prevents its precession near the edge (absence of surface 
charges) [^7], ^0|. Finally, the 0-field can be written 
in the form 



In the case of electrical current influence the circular vor- 
tex dynamics can be excited by the force B7L B8| : 



rg[C-Z] +arg[C-Z J ] -argZ + cJ- (5b) 



Here Z 1 = ZL 2 /R 2 is the image vortex coordinate; the 
image vortex is added to satisfy the Dirichlet boundary 
conditions. 

Under the influence of an external force the vortex 
starts to move. If such forces are weak, the vortex be- 
haves like a particle during its evolution. Such rigid 
vortex dynamics can be well-described using the Thiele 
approach |50| pl] [. Following this approach we use the 
traveling wave ansatz m(r,r) = m(r — R(t)), where 
the function m on the right-hand side describes the vor- 
tex shape (|^). This results in a force balance equation 

F9 + F d + F = 0. The first term F» = -2Q R x z 

is a gyroscopical force [Ell [52|, which acts on a mov- 
ing vortex in the same way as the magnetic field in- 
fluences a charged particle by a Lorentz force. The 
second term is a dissipative force F d = —r)R with 
j] w (a/2) In L, and the last term is an external force 
F = -(27T)- 1 V R S = F ms + F cxt , where the total energy 
is described by Eq. (|J). The magnetostatic part of this 
force is caused mainly by magnetostatic energy of volume 
charges. For small displacements R <C L from the disk 
origin and small aspect ratios (e <C 1), it can be found 
analytically that F ms = -fl G R Q with fl G ~ 10e/97r 
IH . For finite e one can use the analysis in Ref. |33|, |I| , 
which results in £l G (e). Finally the force balance condi- 
tion takes the form of a Thiele equation 
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Without dissipation and external force this equation de- 
scribes a gyration of the vortex around the disk origin 
with a frequency £l G . The dissipation damps the possi- 
ble circular motion of the vortex. However, an external 
force can excite a non-decaying vortex motion. For ex- 
ample, it is known that the rotating magnetic field can 
excite a circular rotation of the vortex in a Heisenberg 
magnet [|o[ |d| and a nanomagnet |[f| due to the force 


F CKt = -^{Rcos(<S>-ljt) + [R x £]sin($-^r)}. (7a) 
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Note that the vortex motion can be excited by any small 
magnetic field (b, lo) p5[ ; howeber, to excite the vor- 
tex motion by a spin current, its intensity should exceed 
threshold value 



NUMERICAL STUDIES 

In order to check the analytical predictions about the 
driven vortex dynamics, we have performed numerical 
simulations. In the case of a rotating magnetic field we 
used the micromagnctic simulator for Landau-Lifshitz- 
Gilbert equations Q as described in Ref. Q . To study 
the vortex dynamics under the action of a spin-polarized 
electrical current, we used numerical simulations of the 
discrete spin-lattice Eq. (^) with the spin torque given 
by Eq. ( ^b| ) as described in Rcfs. [^2|, In both cases, 
as initial condition we use the vortex with positive po- 
larity (p = +1), centered at the disk origin. To identify 
precisely the vortex position we use, similar to Ref. J60| , 
the crossection of isosurfaccs m x = and m y = 0. 

Let us start with the case of a spin current. If we 
apply the current, whose polarization is parallel to the 
vortex polarity (jap > 0), the vortex does not quit the 
disk center, which is a stable point. However, if the spin 
current has the opposite direction of the spin polarization 
(jap < 0, p = +1, a = +1 and j < in our case) the 
vortex motion can be excited when the current intensity 
is above a threshold value |)6], |57], Q as a result of the 
balance between the pumping and damping. Following 
the spiral trajectory, the vortex finally reaches a circular 
limit cycle. The rotation sense of the spiral is determined 
by the gyroforcc, i.e. by the topological charge Q, which 
is a clockwise rotation in our case. At some value, the 
radius becomes comparable with the system size L and 
the vortex dynamics becomes more complicated and can 
not be described by the Thiele equation (^|). Finally it 
results in a switching of the vortex core, see below. 

A qualitatively similar picture is found in the case of a 
rotating field with angular frequency directed opposite to 
the vortex polarity: up < 0. There exist different regimes 
in the vortex dynamics [55| . In the range of frequencies 
close to the frequency of the orbital vortex motion (f2 ~ 1 
GHz) the vortex demonstrates a finite motion in a region 
near the disk center along a quite complicated trajectory 
and does not switch its polarity. These cycloidal vortex 
oscillations are similar to those in Ref. [Q, and corre- 
spond to the excitation of higher magnon modes during 
the motion. If lo 3> ^, then for weak fields the vortex 
motion can be considered as a sum of two constituents: 
(i) the gyroscopic orbital motion as without field, and 
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(ii) cycloidal oscillations caused by the field influence. 
For the case ujp < the direction of the cycloidal oscil- 
lations coincides with the direction of the field rotation, 
while the direction of the gyroscopic orbital motion is 
opposite to it. For stronger fields the vortex motion be- 
comes more complicated and its average motion can be 
directed even opposite to the gyroscopical motion (this 
situation is shown in Fig. [I]). The irreversible switching 
of the vortex polarity can be excited in a specific range of 
parameters (u>, b), with typical frequencies about 10 GHz 
and intensities about 20 mT p5l . The mechanism of the 
switching is discussed below. 



VORTEX POLARITY SWITCHING 

The mechanism of the vortex switching is of general 
nature; it is essentially the same in all systems where 
the switching was observed (2^, |5f| [5(| [n| 
|64| . Under the action of pumping, the original vortex (V, 
qv = l,Pv = lj Qv = — l/2), situated in the disk center, 
moves along a spiral trajectory in the case of a current or 
along a more complicated trajectory (as described in the 
previous section) in the case of the rotating field. During 
its motion, the vortex excites a number of spin waves JS5[ . 
Mainly magnon modes of two kinds are excited in this 
system: symmetrical and azimuthal ones |}6], |3^| . Due to 
the continuous pumping the system goes to the nonlinear 
regime: the amplitude of the azimuthal mode increases 
and there appears an out-of-plane dip nearby the vortex, 
sec Fig. [l](a). When the amplitude of the out-of-plane 
dip reaches its minimum [m z = — 1, Fig. 0(b)] a pair of 
a new vortex (NV, q^v = —^-iPnv = — lj Qnv = — 1/2) 
and antivortex (AV, qAv = 1,Pav = — 1,Qav = 1/2) is 
created. These three objects move following complicated 
trajectories which result from Thiele-like equations. The 
directions of the motion of the partners are determined 
by the competition between the gyroscopical motion and 
the external forces, which are magnetostatic force F ms , 
pumping force F ext , see Eqs. (Q), and the interactions 
between the vortices F- nt . 

The magnetostatic force for the three body system can 
be calculated, using the image vortex approach, which 
corresponds to fixed (Dirichlet) boundary conditions. For 
the vortex state nanodisk such boundary conditions re- 
sults from the magnetostatic interaction, which is local- 
ized near the disk edge |^|. The same statement is also 
valid for the three vortex (V-AV-NV) state, which is con- 
firmed also by our numerical simulations. The 0-field can 
be presented by the three vortex ansatz 

3 

cj> = ^ % {arg[C-^]+arg[C-2/] -arg^J+C^. (8) 

i=l 

The force coming from the volume magnetostatic charge 
density A = — V • m can be calculated in the same way 



as for a single vortex |57| 
-~Vt G q i Y J - ] q]Rj, where & 



|6q ], which results in F ms ss 
= ±1 is the vorticity of i-th 



vortex (antivortex). 

The interaction force F lnt 
coulomb force i* 1 ™* = qiqj- 



between vortices is a 2D 



The internal gyroforces of the dip impart the initial 
velocities to the NV and AV, which are perpendicular to 
the initial dip velocity and have opposite directions. In 
this way the AV gains the velocity component directed 
to the center, where the V is moving. This new-born pair 
is a topologically trivial Q — pair, which undergoes a 
Kelvin motion |67|]. This Kelvin pair collides with the 
original vortex. If the pair was born far from V then the 
scattering process is semi-elastic; the new-born pair is 
not destroyed by the collision. In the exchange approach 
this pair survives, but it is scattered by some angle ]68| . 
Due to additional forces (pumping, damping and magne- 
tostatic interaction), the real motion is more complicated 
and the Kevin pair can finally annihilate by itself, see 
Fig. 0(c). 

Another collision mechanism takes place when the pair 
is born closer to the original vortex. Then the AV can 
be captured by the V and an annihilation with the orig- 
inal vortex happens, see Fig. |l|(d). The collision process 
is essentially inelastic. During the collision, the origi- 
nal vortex and the antivortex form a topological nontriv- 
ial pair (Q = —1), which performs a rotational motion 
around some guiding center (6^, ^9|. This rotating vor- 
tex dipole forms a localized skyrmion (Bclavin-Polyakov 
soliton ffo[), which is stable in the continuum system. 
In the discrete lattice system the radius of this soliton, 
i.e. the distance between vortex and antivortex, rapidly 
decreases almost without energy loss. When the soliton 
radius is about one lattice constant, the pair undergoes 
the topologically forbidden annihilation 55, S9|, which is 
accompanied by strong spin-wave radiation, because the 
topological properties of the system change fnj . 

The three-body problem can be analytically described 
in a rigid vortex approach, based on Thiclc-like equations 

(I): 
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i = 1,2,3. 



Here Ff xt is the driving force (Q) , caused by the field ( |7a| ) 
or spin-current ( f7b| ) pumping. This three-body process 
in exchange interaction approximation is studied in detail 
by Komineas and Papanicolaou p3] . 

The set of Eqs. (||) describes the main features of the 
observed three-body dynamics. During the evolution, 
the original vortex and the antivortex create a rotating 
dipole, in agreement with Refs. |6^|, see Fig. |^. The 
distance in this vortex dipole rapidly tends to zero. The 
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(a) t=110 ps: the rotating field causes the vortex motion, as a result a negative dip is forming. 
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(b) t=320 ps: a vortex-antivortex pair has been born in region 1. 




(c) t=420 ps: the vortex pair has annihilated in region 2, after some time a new vortex-antivortex pair is born in region 3. 
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(d) t=495 ps: the new antivortex has annihilated with the old vortex in region 4 and only the new vortex of opposite 

polarity remains. 

FIG. 1: The vortex switching process under rotating field influence. Left column - vortex trajectories. Red line - initial 
vortex core trajectory, green line - new vortex, blue line - antivortex. Right column - in-plane magnetization distribution, the 
out-of-plane component is denoted by levels of gray color, dashed blue and solid green lines correspond to isosurfaces ra, = 
and m y — 0, respectively. The out-of-plane component m z is shown in the middle column. The data were obtained from 
micromagnetic simulation for a permalloy disk (132nm diam., 20nm thickness) dynamics under rotating field (B = 0.07 T, 
uj = 10 GHz) influence. 
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FTC! ?■ Ttirep-hody dynamics from the solution of Eg (|p]), 
for £ system with LI I = 20. a = 0.01. & _ 0.02. uj = 1. and 
JIg j= 0.0287. The original vortex (red line) was positioned 



at (1,0)£, the antivortex (blue) at (2, 
vortex (green) at (2.5, 0)1. 



and the newborn 



new-born vortex moves on a clockwise spiral to the ori- 
gin. 
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SUMMARY 



To summarize, we have studied the switching of the 
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